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Abstract 



«4H _ 

O I In this note I present a simplified version of tlie recent computation (Mezard and Parisi 

1998, 1999) of the properties of glasses in the low temperature phase in the framework of the 
^2 ' replica theory, using an extension of the tools used in liquid theory. I will only consider here the 

case of the internal energy at T = 0, which can be studied in a simple way without introducing 
' replicas. 

a ■ 
o 

1 Introduction 

^ , In this note I will present the recently developed theory, based on replicas, (Mezard and Parisi 
I 1998,1999) which aims to derive from first principles the basic thermodynamic properties of glasses 
^ using a modified liquid theory. In order to simplify the presentation I am studying only one point: 
•/"J ■ the computation of the ground state energy, i.e. the internal energy at zero temperature {Em)- To 
this end I am going to use an explicit approach, in which the physical assumptions are evident. In 



G^ , this approach no replicas are introduced, but the same results can be obtained using the replica 
method. This note can be considered a pedagogical introduction to the more sophisticated approach 
of Mezard and Parisi (1998, 1999). 
r ■ Beyond technicalities the basic assumption is that a glass is very near to a frozen liquid. This 
quite old statement can be rephrased as a liquid is very near to a heated glass. In other words the 



O . configurations of a glass at low temperature are not far from those of a liquid. Therefore if we use 
'P. ' some smart method to explore the phase space in the liquid phase we can find the properties in the 

low temperature glassy phase. 
^ This strategy has been put into action using the replica theory. In this way one finds, as output 
' of an explicit computation, that the glass transition is characterized by the vanishing of the configu- 
rational entropy (the so called complexity) and in the replica language the low temperature phase is 
described by one step replica symmetry breaking (Mezard et al. 1987, Parisi 1992, Kirkpatrick and 
Thirumalai 1995) with a non vanishing non-ergodic parameter at the phase transition point. The 
computations can be done at all temperatures. Here we are going to reproduce the results of the 
replica approach at T = 0, where only energetic considerations (not entropic ones) are relevant. 



2 Basic assumptions 

The strategy of the computation is rather simple. Given a system of particles with an Hamiltonian 
if(x), we would like to know the minimum value of H. This task in non trivial when ^ cxd if the 
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ground state is not translationally invariant. The minimum can be characterized as the solution of 
the stationary equations 

having the minimum energy density E = H/N . 

In order to put the problem in a wider perspective it is convenient to consider the function 

N 

u{E)=Y.5{E-E^)= d^\dei{n{-K)\5{NE - H{-K))\{5{Fi{:s)), (2) 

Q ■' i = l 

where the sum runs over all the solutions of eq (|Tp and the Hessian matrix 7i is an iV x iV matrix 
defined as 

Ti-i^k = • (3) 

There are many arguments which imply that in the limit N going to infinity we have (in distri- 
bution sense) 

iy{E) ^ exp{Na{E)), (4) 

where the function (y{E) is positive for E > Em, and vanishes a.t E = Em'- ^{E) = for £" < Em- 
Therefore our goal is reached if we evaluate the function (j{E). 

Unfortunately the computation of may be too complicated. In many cases it is convenient 
to define the function 



fi{E) = Y,S{E- Ea) sign(det(H)(x))) = / rfx det{n{^)d{NE - i7(x,)) H (5(F,(x)). (5) 

The advantage of this second definition is that the final integrand is a smooth function of x and 
this may simplify the computation. 

In the same way we define in the large limit the function 

fi{E) ^ exp{NE{E)). (6) 

Of course in the region where the Hessian is positive the two definitions coincide. 

In the region of E near to Em we expect that most of the stationary points are minima and there- 
fore the function ct{E) and ^{E) are similar in this region (Cavagna et al. 1998). As a consequence 
both functions can be used for the computation of the ground state. In the rest of the paper we are 
going to compute S(-E), however within the approximation we will do in this note the two functions 
coincide. 

The evaluation of ^{E) may be simplified (Monasson 1995) if we introduce the generalized par- 
tition function 

N 

C(7) = J dEexp{-N{^E-J:{E))) = j d^ det{H{^)) exp{-^H{^))Yl6{F,{^)). (7) 

i=l 

In the large N limit we have that 



2 



C(7) = cxp(-iV<l>(7)), 
$(7)=7£;(7)-S(£;(7)). (8) 
a$(7) 



^7 



The previous equations are the equivalent of the usual relation among the free energy and the 
entropy. The quantity E is the equivalent of the entropy and is called the complexity or the config- 
urational entropy. 

We now face the problem of computing the quantity $(7). We are clearly interested in the 
region of relatively large 7, which corresponds to small energies and consequently small S. The 
computation of $(7) is rather involved. In order to simplify it we can make three assumptions on the 
configurations which dominate the integral at large 7, that are similar to the configurations which 
dominate the usual partition function at large /3. 

• The integrand is dominated by those configurations in which all the eigenvalues are positive, 
i.e by configurations near to a minimum, not to a saddle point. 

• The dominant integration region splits into many, approximately disjoint, regions (TZa) and in 
each of these regions the Hamiltonian can be well approximated by a quadratic one. 

• The variations of the determinant of the Hessian from one region to an other region can be 
neglected. Therefore the determina nt of the Hessian is nearly constant in all the regions. 

3 How to use the assumptions 

We can now proceed. We use the following representation of the 6 function: 

S{Fk) = I dXkeMi>^kFk). (9) 

We thus find 

C(7) = / dKdX det(?i(x)) exp ^-^H(x) + i A^F^ j . (10) 

If we now shift the integration over the x variables and we introduce the variables i/k — Xk+i'y''^Xk 
we obtain 

C(7) = / dydXdet{n{y))exp ||_^i^(y) _ |^ A,A,H(y),,, + O(i/^^oj ■ (H) 

We can now neglect the anharmonic terms proportional to the third derivative of H, (i.e. the 
terms of order H^^^) to perform the integral over the A. We arrive at the simpler expression 



C(7) = / o?yexp(-7i/(y) + $^(7,y)), (12) 



where 
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exp(-<l>^(7,y)) = I dXexp ^--fj^ \k\iniy)k,}j 



f— j detin)-'/^ = exp(-l/2Tr(ln(7i:))). (13) 

The formula (12) can be derived directly from the previous assumptions. Indeed if we start from 
eq. ([T^) and neglect the anharmonic terms in each region, we find that the determinant is constant 
so it can be taken out of the integration and we obtain: 

^exp(-<|.^(7)) / dyexpi-^Hiy)) = Y^expi-^E)^. (14) 

Indeed the integral in each of the regions TZa in which the integrand is approximately harmonic 
gives approximately exp{—'yEo, + $^^(7)) (Angelani et al. 1998). 

Although the formula (|12D can be derived directly, the previous detailed proof is useful in order 
to understand which are the possible steps that must be followed in order to improve its validity if 
some of the assumptions are removed. 



Also after these simplifications the evaluation of eq. (12) is not trivial because the determinant 
changes from point to point. Assuming that it takes the same value near all the low lying minima, 
we get 

C(7) = ^(7) exp ($^(7)) = exp (-7^(7) + $^(7)) , (15) 

where $^(7) is the average value of the harmonic free energy (multiplied by 7) for a configuration at 
temperature I/7 and Z^j) is the usual partition function at /3 = 7. We need to know only the free 
energy (-^(7)) and the spectral density of the INN (instantaneous normal modes, i.e. the eigenvalues 
of the Hessian) in the liquid phase as a function of the temperature in order to find C(7)- 

Let us use the previous formulae for $^(7) and let us introduce the corresponding harmonic 
entropy: 

^^(7) = -$^(7) + ^^~ ln(7) + A(7), (16) 

where A(7) depends on 7 only via the eigenvalues of the Hessian in the liquid at a temperature 7"^. 
After a simple computation we find that the complexity is given by 

E(7) = S^(7)-^^(7)+7^, (17) 

where S^{'y) is the entropy in the liquid. 

In this way we get a simple expression for the complexity in terms of quantities that are defined 
in the liquid. We need to known then in the high temperature phase down to a temperature I/7*, 
where ^(7*) = 0. 

The zero temperature configurations of the glass will not be the same as those of the liquid 
at finite temperature. In the harmonic approximation we have to find the minima by assuming a 
quadratical energy and estimating the position of the particles in a glass using the formula 

y.=Xi + 6i, (18) 

where Xi and yi are the position of the particles in the minimum and in the liquid respectively and 

^ n,,A = F,. (19) 

k=l,N 



The analytic solution of the equation (|T9D is not so simple, and it will not be discussed here. 
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4 How to implement previous ideas 



It is always possible to evaluate eq. (0) using numerical simulations in the liquid phase (Coluzzi 
et al 1999a and 1999b). This is by far the simplest method. I now present a possible analytic 
implementation of this strategy (Mezard and Parisi 1998, 1999). 

The computation of the free energy in the liquid phase does not present any serious difficulty. 
There are many methods which lead to integral equations for the correlation functions (e.g. the 
hypernetted chain approximation, the mean spherical approximation, the Percus Yevick approxima- 
tion). We have used the hypernetted chain approximation, however the other approximations may 
work better. In any case a clean computation of the properties of the liquid is standard (in the case 
of two body central interaction). 

The computation of the spectrum is much more difficult. Most of the computations are done at 
low density and in spite of this simplification the results obtained are approximated (for example 
the localized tail at small and large eigenvalues is lost). The only exact computation at low density, 
which corresponds to the first order in the virial expansion, has just been done recently (Cavagna et 
al 1999). Lacking for the moment a systematic way to extend this computation to the high density 
region, we are following a quite simple method which becomes exact in some limit. A closed formula 
can be obtained if we assume that the range of the Hessian is much larger than the typical interatomic 
distance and we use the superposition approximation 

g{xiX2X3) = g{xi - X2)g{x2 - x^)g{x^ - xi) (20) 

in the cases where it is needed (Mezard and Parisi 1998, Mezard et al. 1999). 

In order to compute the eigenvalues of the Hessian we compute analytically the moments of the 
eigenvalue distribution 

Mk = Ti^rt''). (21) 

Each moment can be written as the appropriate integral over the correlations functions. In the 
long range approximation some terms are dominant over the others. We keep only those terms and 
we introduce the superposition approximation into the appropriate places. We use the large range 
expansion to select the diagrams. The result we obtain for the spectrum is just a first approximation, 
the long range expansion is a nice way to organize the diagrams and we plan to improve it in the 
next future. An unified approach, which interpolates in a reasonable way the expressions in the low 
density with those in the high density limit would be very useful. 

After some computations (which we do not report here) we find a spectral density p(e) that has 
support only in the region of positive eigenvalues e (real frequencies u = \J{e)) and goes to zero as 
uo as expected. Some extra work is needed to improve the approach and having a realistic S{k,uj) at 
T = 0. 

We have applied this formalism to soft spheres with repulsive interaction (Mezard and Parisi 1998, 
1999), to binary mixtures of soft spheres (Coluzzi et al. 1999a) and to binary mixture of particle 
with LJ potential (Coluzzi et al. 1999b). We plan to study in the future undercooled Argon and 
81204. 

Here I present the results at T = for soft spheres with a repulsive x~^'^ potential. In fig. ^we 
see the complexity as function of free energy at three different temperatures (at zero temperature the 
free energy become identical to the internal energy). The function S should be considered reliable at 
low energy, however at high energy the harmonic approximation fails and non linear terms become 
more important. Moreover in the real world increasing the energy the number of negative (maybe 
localized) eigenvalues increases and our approximation becomes no more reliable. 
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Figure 1: The configurational entropy versus the free energy, (from Mezard and Parisi 1999) 
computed within the harmonic resummation, at temperatures T = 0, .05, .1 (from left to right). 

If we extend the same arguments from energy to free energy (following Monasson 1995) we must 
introduce replicas. We find that there is transition at temperature Tk such that: 

• For temperature greater than the number of minima which can be reached by slowly cooling 
the equilibrium configurations increase as exp{NT:{E,T)). The equilibrium configurational 
entropy goes to zero at Tk- 

• For temperature smaller that there is a finite number of of minima which can be reached by 
slowing cooling the equilibrium configurations. The configurational entropy is zero below T^. 

• The transition temperature corresponds to a value jk which satisfies the following relation 
(Mezard and Parisi 1998) 

E^(7) = S^(7) (22) 

We notice that, if we neglect the temperature dependence of A(7), this equation coincides 
with the condition E(7) = 0. The value of coincides with the value of 7* which is used to 
compute the ground state energy. 

This is an explicit realization of a transition driven by an entropy crisis, which is the first simple 
realization in the random energy model (Derrida 1980). The specific heat has a jump downward, 
when we decrease the temperature, which is the opposite of the typical behaviour in transitions 
characterized my the onset of a conventional order parameter (ferromagnets, superconductors, have 
a jump upward when we decrease the temperature). 

5 Conclusions 

We have found a method which is able to use liquid theory method in the glasses phase putting in 
practice the old adagio a glass is a frozen liquid. 

We are able to compute with a reasonable approximation the thermodynamics and with a little 
more effort we can compute the static and the dynamic structure functions. 

I would finally stress that in this theoretical framework we do not only obtain reasonable results 
for conventional experiments, but we have also clear and quantitative predictions for the relations 
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between the response and the correlations in the aging regime (Cughandolo and Kurchan 1993, Franz 
and Mezard 1994, Franz et al. 1998, 1999). An experimental test of these relations, which can be 
done with present technology but with a lot of ingenuity, would have extremely deep consequences 
for theory (this theory would be just killed in the case of a negative result, while the whole approach 
based on replicas would be confirmed in the positive case). These violations have been seen quite 
clearly in numerical simulations (Parisi 1997a and 1997b) and the time is ripe for seeing them in real 
experiments. 

I hope that the experimentalists will consider such a measurement as a challenge and the results 
will be available for the next Andalo meeting. 
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